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Antiferromagnetically ordered state of monoclinic quasi-one-dimensional S = 1/2 Heisenberg antiferro¬ 
magnet CuSb 206 was studied combining torque magnetometry with phenomenological approach to magnetic 
anisotropy. This system is known to have a number of different twins in monoclinic (3 phase which differ in ori¬ 
entation of the two CUO 5 octahedra in unit cell resulting in different orientation of magnetic axes with respect to 
crystal axes for each twin. We performed torque measurements in magnetic fields H < 0.8 T on a sample where 
certain type of twin was shown to be dominant by ESR spectroscopy. The measured data reveal that easy axis is 
the crystallographic b axis for this sample. Phenomenological magnetocrystalline anisotropy energy invariant to 
crystal symmetry operations was used to model spin axis direction in zero and finite magnetic field. Our model 
reproduces the value of the spin flop field H^f = 1.25 T found in literature. Combination of this approach 
with our torque results shows that the spin axis will flop in the direction of the maximal value of measured g 
tensor when magnetic field H > Hgp is applied along easy axis direction. Our analysis of magnetocrystalline 
anisotropy energy predicts two possibilities for the easy axis direction in this system, b or a, connected to differ¬ 
ent crystallographic twins that can be realized in CuSb206. These results offer possibility to reconcile different 
reports of easy axis direction found in literature for this system and also nicely demonstrate how combination 
of torque magnetometry and phenomenological approach to magnetic anisotropy can be used to determine the 
value of the spin flop field and the direction of spin axis in antiferromagnets in both H < Hgp and H > Hgp by 
performing measurements in fields significantly smaller than H^p- 

PACS numbers: 75.25.-j, 75.30.Gw, 75.50.Ee 


I. INTRODUCTION 

Copper oxides represent one of the most studied materials 
in condensed matter physics. While this has its roots in real¬ 
ization of first high-temperature superconductors in cuprates, 
nowadays these materials present an arena for study of many 
interesting and potentially applicable phenomena related to 
magnetism. Magnetism in copper oxides comes from copper 
ion Cu^+ which has unpaired electron carrying spin S = 1/2. 
Magnetic Cu^+ ions are treated as spin-only and the effect of 
their unquenched orbital moment is included in anisotropic 
electron g tensor. Distorted and usually elongated CuOe 
octahedra found in most copper oxides are responsible for 
orbital ordering of the 3(1 orbital of the unpaired electron 
which then influences the dimensionality of the specific 
magnetic lattice defined by the dominant interactions be¬ 
tween spins S = \/2 and usually described by the isotropic 
Heisenberg Hamiltonian iH = J'E.ij^i '^j- In magnetically 
ordered states, however, magnetic moments choose specific 
directions in space which are not anticipated by this isotropic 
model. Weak anisotropic terms must then be added to the 
spin Hamiltonian which often significantly complicates the 
theoretical treatment of the properties of these systems. Exact 
experimental determination of magnetically ordered ground 
state can be helpful in identifying the relevant anisotropic 
contributions in specific system. 
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In this work we study the antiferromagnetically ordered 
state of monoclinic quasi-one-dimensional CuSb206 mo¬ 
tivated by the different reports on easy axis direction and 
ground state magnetic structure found for this system in 
literature. We employ torque magnetometry as a sensitive 
probe of macroscopic magnetic anisotropy. Main advan¬ 
tage of torque magnetometry over magnetic susceptibility 
measurements is that magnetic torque measures in-plane 
magnetic susceptibility anisotropy rather than susceptibility 
along specific axis. This allows not only determination of the 
direction of easy axis in collinear antiferromagnet, but also 
predictions about the directions in which the spins will flop 
when magnetic field along easy axis direction is applied. Our 
analysis of the interplay of crystal structure and magnetic 
anisotropy of different crystallographic twins realized in 
CuSb206 in combination with our torque measurements puts 
forward a possibility to reconcile different reports of the 
ground state magnetic structure of this system. 

CuSb206 crystallizes in monoclinically distorted trirutile 
structure type [1] at room temperature. Below 380 K a second 
order phase transition takes place from high-temperature 
tetragonal a—CuSb206 (space group P42/nmn) to mono¬ 
clinic p—CuSb206 (space group P2\/n) [2, 3]. This phase 
transition lifts the degeneracy of four 02 oxygens which 
exists in the tetragonal phase. Above 380 K dynamic 
Jahn-Teller effect is realized in the tetragonal phase where 
higher vibronic states are occupied leading to fast dynamic 
exchange between the two different possibilities for elonga¬ 
tion resulting in averaged compressed CuOg octahedra [2]. 
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Below the phase transition static Jahn-Teller effect is realized 
through localization of vibronic states in two energy minima 
of the ground-state potential surface, corresponding to CuOe 
octahedra elongated along the Cu-02a or the Cu-02 bond 
[2] . The phase transition from tetragonal to monoclinic phase 
is accompanied by a formation of crystallographic twins 
which is a result of two different choices for elongation of 
two crystallographically equivalent CuOg octahedra in the 
unit cell of CuSb206 [3, 4]. 

Magnetic susceptibility measurements revealed that 
CuSb 206 is quasi-one-dimensional antiferromagnet with 
intrachain exchange interaction J/ks ~ 90 K [1], which 
was unexpected given the almost square planar arrangement 
of magnetic Cu^+ ions in the ab plane. Octahedral oxygen 
environment of the magnetic Cu^+ ions gives rise to dominant 
Cu-O-O-Cu superexchange pathway which forms magnetic 
chains along the {a + b) direction at z = 0 and {a — b) direc¬ 
tion at z = 0.5 [5, 6]. Band structure calculations interpreted 
the origin of one-dimensional magnetism in CuSb 206 as 
being driven by unusual orbital ordering which allows su¬ 
perexchange to be realized through dj^ 2_^2 orbitals and apical 
instead of square planar oxygens [7]. These calculations, 
however, were performed for tetragonal phase of CuSbzOg 
where CuOg octahedra are compressed. In monoclinic phase 
standard orbital ordering of d^ 2 _y 2 orbitals is expected due to 
the above mentioned elongation of octahedra [6]. 

Despite the one-dimensional character of the magnetic 
lattice of CuSb 206 , phase transition to long range antiferro¬ 
magnetic order takes place below TX? w 8.5 K due to weak 
interchain interaction [4, 5, 8, 9]. The reports on magnetic 
structure and magnetic susceptibility anisotropy are, however, 
somewhat controversial. Neutron diffraction measurements 
determined the propagation vector k = (1/2,0,1/2) and 
reduced magnetic moment of ^eff ~ 0.5/jb [5, 8, 9] and 
11 ^ff « 0.4/jb in one study [10]. Magnetic moments within 
the ab plane are oriented antiparallel along a and c directions 
and parallel along b direction. Different orientations of 
magnetic moments with respect to crystal axes are, however, 
reported in literature. Powder neutron diffraction measure¬ 
ments revealed spins to be confined to ab plane, but could 
not distinguish between structure where spins at z = 0 and 
z = 0.5 are collinear and structure where they are canted 
[5]. Single crystal neutron diffraction studies in one case 
proposed magnetic structure with moments canted away from 
the b axis toward the a axis with different canting in z = 0 
and z = 0.5 layers [9]. This choice of canting was based on 
the magnetic susceptibility anisotropy reported in the same 
paper [9]. Two other single crystal neutron diffraction studies 
showed that collinear antiferromagnetic order is realized with 
spins aligned in the b axis direction, although some weak 
reflections could not be indexed [8], or were disregarded as 
signal from crystal twins [10]. 

Below Tn different magnetic susceptibility anisotropy 
reports are found in literature. Large decrease of suscep¬ 
tibility characteristically observed when magnetic field is 
applied along the easy axis has been reported for b axis 
direction [4, 11], a axis direction [3], and even for both 
directions [9, 12, 13]. Recent report on single crystal mea¬ 


surements of magnetic susceptibility showed that significant 
decrease of susceptibility below TV and spin flop transition 
at Hsf = 1.25 T is observed when magnetic field is applied 
along both the a and b axis, suggesting that both a and b are 
the easy axes in CuSb206 [12]. The magnetically ordered 
ground state of CuSb 206 is thus still not completely resolved. 

In this work we attempt to shed some light on magnetically 
ordered ground state of CuSbzOe using torque magnetometry 
in magnetic field smaller than the spin-flop field. We combine 
our torque measurements with phenomenological approach to 
magnetic anisotropy where we model the magnetocrystalline 
anisotropy using crystal and magnetic symmetry restrictions 
to describe the anisotropic magnetic response of this system 
in finite magnetic held. This approach allows us to determine 
the orientation of spin axis in both zero held and in // > Hsf- 
In Sec. II we describe the characterization of our sample by 
X-ray diffraction, magnetic susceptibility measurements and 
electron spin resonance (ESR) spectroscopy. In Sec. Ill A the 
results of torque magnetometry are presented. Phenomeno¬ 
logical approach to magnetic anisotropy of CuSb 206 is used 
for the description of measured torque curves in Sec. IIIB. 
The results are discussed and the conclusion given in Secs. 
IV and V, respectively. 


II. SAMPLE CHARACTERIZATION 

Single crystals of CuSb206 were grown using endothermic 
chemical vapor transport reactions with different transport 
agents. Firstly, polycrystalline CuSb206 powder was synthe¬ 
sized by solid state reaction described in detail elsewhere [3]. 
Secondly, CuSb206 powders and HCl or TeCU were mixed 
and sealed in a 24 mm diameter and 20 cm long evacuated 
quartz tube. The tube was inserted in a two-zone horizontal 
tube furnace with a source-zone temperature of 900°C and 
a growth-zone temperature of 800°C for four weeks. In 
this way, millimeter-sized single crystals of CuSbzOe were 
obtained. No impurity phases were detected using X-ray 
powder diffraction. The crystals are yellow in color, as 
reported in literature [3]. 

One large crystalline sample was ground into powder and 
used for measurements of average magnetic susceptibility. 
Small thin crystal with mass m = (91 ±2) pg cut from 
larger crystal was used for torque measurements. On the 
surface of the sample thin stripes due to twinning were visible 
[3]. We did not attempt to detwin the sample. This sample 
was characterized by X-ray diffraction and ESR spectroscopy. 


A. X-ray diffraction 

The orientation of the crystal axes with respect to the 
sample morphology was determined using X-ray diffraction. 
The crystal was glued onto a glass fiber and the data were 
collected on an Oxford Diffraction Xcalibur 3 CCD diffrac¬ 
tometer with graphite-monochromated Mo Ka radiation 
(A. = 0.71073 A) at room temperature. The orientation of 
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FIG. 1. Magnetic susceptibility of powder CuSb206. Solid line rep¬ 
resents fit to magnetic susceptibility of spin S = 1/2 ID Heisenberg 
antiferromagnet [15, 16]. 


FIG. 2. X-band ESR spectra (right) and g factor anisotropy (left) 
measured in the be plane at T = 80 K. Solid lines in figure on the right 
represent results of fitting of spectra to four individual Lorentzian 
lineshapes. Error bars for g factors obtained from the fit are smaller 
than the size of symbols. 


the crystal and the unit cell parameters were obtained by the 
CrysAlis software package by Agilent Technologies. Data 
confirmed the monoclinic space group P2\jn with the follow¬ 
ing unit cell parameters: a = 4.6299(3) A, b = 4.6343(3) A, 
c = 9.2911(6) A, p = 90.965(7)°, in good agreement with 
previously published data [1]. Presence of multiple twins, 
reported by other authors [2-4, 12-14], was also observed in 
our data. 


B. Magnetic susceptibility 

Magnetic susceptibility of powder, {%), was measured us¬ 
ing Faraday apparatus in the temperature range from 2 K to 
330 K in magnetic field H = 1 kOe. Temperature depen¬ 
dence of magnetic susceptibility (%) of CuSb 206 measured 
on powder is shown in Fig. 1. The measured data can be 
described by the susceptibility of one-dimensional (ID) spin 
5=1/2 Heisenberg antiferromagnet with J/ks = 94 K and 
{g) = 2.11 [15, 16], in agreement with results reported in lit¬ 
erature [1, 3,4, 8,9, 12]. At Tn ~ 8.5 K akinkin susceptibility 
marks the phase transition to long range magnetic order. Be¬ 
low 7// (x) decreases as temperature decreases in a manner 
typical for collinear antiferromagnet [17]. To determine the 
direction of the easy axis we employ torque magnetometry in 
Section III A. 


C. ESR spectroscopy 

X-band ESR measurements were performed at T = 80 K 
using Bruker Elexsys 580 ET/CW X-band spectrometer 
equipped with a standard Oxford Instruments model DTC2 
temperature controller. The microwave frequency was ^ 
9.7 GHz with the magnetic field modulation amplitude of 5 G 
at 100 kHz. ESR spectra were recorded at 5° steps and the 
rotation was controlled by a home-built goniometer with the 
accuracy of 1°. The uncertainty of w 5° was related to the op¬ 


timal deposition of the crystal on the quartz sample holder. 

We performed ESR spectroscopy measurements at T = 
80 K to determine the number and type of twins in our sin¬ 
gle crystal sample which is used for torque magnetometry 
study described in Sec. Ill A. Previous detailed ESR analy¬ 
sis of CuSb 206 has shown that there are four possible twin 
domains in the monoclinic p phase which appear simultane¬ 
ously during the structural phase transition [2, 4]. 

Spectra measured at several different directions in the be 
plane of our sample are shown in Eig. 2. Eour ESR peaks 
of different intensity could be detected while for H\\b only 
one single Lorentzian line is observed. Two central lines, 
one of which was the most intensive line in the spectra (twin 
1) and the other was weak (twin 2), are located at approxi¬ 
mately the same resonant fields and therefore it was difficult 
to distinguish them. Each of the four lines can be explained 
by the anisotropic Zeeman interaction due to the anisotropic 
g tensors of the two magnetically inequivalent Cu^+ ions in 
the unit cell which merge to one resonance due to strong ex¬ 
change coupling [4] . To obtain angular dependence of g factor 
for each line, when possible, we fitted our spectra to four in¬ 
dividual Lorentzian lines parametrized by its resonance field 
Hres g, linewidth AH and intensity. The results of fit are 
shown as solid lines describing the spectra in Eig. 2. Due 
to low resolution of our X-band spectra, the contributions of 
twins 1 and 2 were difficult to distinguish at almost all angles, 
and contributions of all four twins were difficult to distinguish 
around b and c axes where spectra of different twins merge. 
Intensities resulting from the fit suggest that twins 3 and 4 
make up for less than 10% of our sample. 

Appearance of the observed four different twins in p- 
CuSb206 has its origin in a change occurring in CuOe octahe- 
dra as they go through phase transition from high-temperature 
tetragonal to the low-temperature monoclinic phase. There 
are two differently oriented but crystallographically equiva¬ 
lent CuOg octahedra in the unit cell of CuSb206. In the tetrag¬ 
onal a phase CuOe octahedra are compressed with the direc¬ 
tion of compression being along the Cu-Ol bonds, and four 
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^\\ea ^ ^SF 

m\ 

m2 

m 3 



Cu-O in Octl 

Cu-O in Oct2 




axis 

I 

P12y/nl 

Cu-02a 

Cu-02a 

(Y,0,r|) 

b 

(-il.O.Y) 

b 

m2, 

II 

Cu-02 

Cu -02 

(Y,0,-ti) 

b 

(t1,0,y) 

b 

m2 

III 

P2i/nIU 

Cu-02a 

Cu -02 

(0,-Y,ri) 

d 

(0,ri,Y) 

d 

m2 

IV 

Cu-02 

Cu-02a 

(0, -Y, -r|) 

d 

(0,-r|,Y) 

d 

m2 


“ This choice of space group and consequently direction of easy axis is obtained if common crystal axes are assumed for all twins, as explained in the text. 


TABLE I. Elongated Cu-O bonds and directions of macroscopic magnetic axes my, m 2 and m 3 with respect to the crystal axes for each twin. 
my, m 2 and m 3 are unit vectors. Octl (Oct2) is lower orange (upper blue) octahedron in Fig. 3. Direction m 3 is determined as the direction 
of the maximal value of the g tensor for each twin. 7 = cos(<(m 3 ,c)) and r| = y/l — Easy axis direction for each twin is given in column 
before last. The direction of spin axis when field > Hsf is applied along the easy axis direction is m 3 in our model, where m 3 axis for 

each twin is given in the second column before last. 



FIG. 3. Two magnetically inequivalent CuOe octahedra in twin I of 
P—CuSb206 as viewed in the ac and the be plane. Lower (orange) is 
represented as Octl in table I, and upper (blue) as Oct2. Crystal axes 
directions are shown in black, and macroscopic magnetic axes my, 
m 2 and m 2 for twin I are shown in red, green and blue, respectively. 


equivalent Cu-02 bonds in the equatorial plane. In the mon¬ 
oclinic p phase CuOe octahedra become elongated along two 
of the Cu-02 bonds which now become Cu-02a bonds, while 
the other two Cu-02 bonds shrink to form an equatorial plane 
with Cu-Ol. Direction of compression in a phase and elonga¬ 
tion in p phase is thus not the same. Four different twins can 
emerge from two different choices for elongation for each of 
the two magnetically inequivalent CuOe, as already discussed 
in Ref. 2 and evidenced in ESR measurements [4]. Using ter¬ 
minology of Ref. 4, we label those twins in roman letters from 
I to TV. In Fig. 3 we show the two CuOe octahedra for twin I, 
Octl and Oct2, in orange and blue. In table I we summarize 
the elongated Cu-O bonds for four possible crystallographic 
twins. 

From the above presented analysis and the results of Ref. 
4 it follows that lines of twins I and II should be indistin¬ 
guishable in the be plane. The splitting we observe for our 
twins 1 and 2 could result from the misorientation of the sam¬ 
ple. The values of g factor we obtain for twin 1 in the be 
plane are gb =2.21 and gc = 2.21. For twin 2 we observe 
gc = 2.24. This compares rather well to the data from litera¬ 
ture, gb = 2.21, gc = 2.23 given in Ref. 4. Our twins 3 and 4 


correspond to twins III and IV in Ref. 4, respectively. 

Local magnetic axes, x, y, and z of each CuOg octahedron 
are defined by the CuOq plaquette of four nearest oxygens 
[2, 4]. Local X axis is in the direction of Cu-Ol, y is in the 
plaquette, nearly in the Cu-02 direction, and z is perpendic¬ 
ular to the plaquette [2, 4]. z axis is thus « 12° away from 
the apical Cu-02a bond of the octahedron. This choice of lo¬ 
cal magnetic axes was corroborated by the measured g factor 
anisotropy for each twin [4]. The resultant macroscopic mag¬ 
netic axes of the sample, my, m 2 and m 3 , are thus not in the 
direction of the crystal axes. Direction of macroscopic mag¬ 
netic axes in Cu^+ 5=1/2 systems can be determined ex¬ 
perimentally as the axes in which the tensor and magnetic 
susceptibility tensor x are diagonal. 

Using the above given local magnetic axes we obtain for 
twin I that macroscopic magnetic axis m 2 is in the direction 
of the b axis, m 2 = b, and my and m 3 are in the ae plane. 
The angle between axes my and c in twins I and II amounts 
to <{my,c) « 49°, and <{mj,c) ~ 41°, but the orientation of 
these axes is different in those twins, as evidenced from LSR 
measurements [4]. In twins III and IV b axis is no longer mag¬ 
netic axis, as was observed in previously reported LSR mea¬ 
surements [4], and as can be seen from our results for twins 3 
and 4 shown in Lig. 2. This has consequences for the crystal 
symmetry, as we will now argue. In twins I and II b axis is 
one of the magnetic axes and the other two magnetic axes are 
in the ae plane and are rotated from a and e crystal axes, as 
evidenced from LSR measurements [4]. According to Neu¬ 
mann’s principle any physical property of a crystal must be 
invariant to the same point group symmetry operations as the 
crystal structure [18]. Lor CuSb 206 this means that measured 

tensor and susceptibility tensor must be invariant to 2 /m 
symmetry operations [19] with respect to b axis for conven¬ 
tional choice of b axis as unique axis in P2y/n space group, 
which can also be written as P12i/nl. Lor twins I and II that 
is indeed the case. Lor twins III and IV b axis is no longer 
magnetic axis, as can be seen from Lig. 2 and also in Ref. 4. 
Magnetic axes for those twins are axis m 2 = d, and axes my 
and m 3 which correspond to directions of minimal and maxi¬ 
mal value of g in the be plane. Magnetic properties such as g^ 
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tensor and magnetic susceptibility tensor are now invariant to 
2/m symmetry elements with respect to the a axis. According 
to Neumann’s principle this means that unique axis is a axis 
for twins III and IV and their space group should be labeled 
as /’2i /n 11. We distinguish these two symmetry variations of 
P2\/n in table I. This will be helpful in construction of mag¬ 
netocrystalline anisotropy energy for each of the twins in Sec. 
IIIB. 

Another possibility of interpreting g factor data for twins 
III and IV would be by assuming that these twins have axes a 
and b rotated by 90° with respect to those of twins I and 11. We 
are not able to either support or dispute this possibility from 
our structural measurements and since it was not suggested 
by previous structural measurements [2-4] we continue by as¬ 
suming common crystal axes for all twins in our sample. 

III. RESULTS 

A. Torque magnetometry in the AFM state of CuSb206 

Magnetic torque was measured with home-built torque ap¬ 
paratus which uses torsion of thin quartz fiber for torque mea¬ 
surement. The sample holder is made of ultra pure quartz and 
has an absolute resolution of 10^“* dyn cm. The uncertainty of 
Ri 5° was related to the optimal deposition of the crystal on the 
quartz sample holder. This uncertainty was somewhat larger 
for the a*b plane because of sample’s morphology. 

Magnetic torque r acting on a sample of volume V with 
magnetization M in magnetic field H is given by expression 

x = yMxH (1) 

In our experiment magnetic field is rotated in an xy plane and 
only Tj component of torque perpendicular to that plane is 
measured. In case of linear response of induced magnetization 
to applied magnetic field, M = x • H, where x is magnetic 
susceptibility tensor, we obtain for the measured component 
of torque 

''' = - %>-) sin 20. (2) 

2M,nol 

In expression (2) m is the mass of the sample, M^oi is 
molar mass, X-x nnd Xv are susceptibility tensor components 
along directions x and y expressed in emu/mol and 9 is the 
angle magnetic field makes with the x axis. Directions x and 
y represent directions along which the torque is zero, i.e. 
magnetization is parallel to applied magnetic field. These 
directions are the directions of magnetic axes in the xy plane. 
From expression (2) we see that measured component of 
torque is proportional to magnetic susceptibility anisotropy, 
^Xxy = Xx — Xy, in the plane of rotation of magnetic field. This 
result is applicable to paramagnets and also collinear antifer- 
romagnets when applied magnetic field is much smaller than 
the spin flop field Hsf- Measurements in the paramagnetic 
state (not shown here) were performed on a larger sample. 
Angular dependence obeys Eq. (2) and magnetic susceptibil¬ 
ity anisotropy displays the same temperature dependence as 



FIG. 4. Temperature dependence of magnetic susceptibility 
anisotropy measured in ac, a*b and be planes. 



FIG. 5. Angular dependence of torque x multiplied by 2M„jgi / (mH^) 
[see Eq. (2)] measured in if = 8 kOe in the ac plane. Macroscopic 
magnetic axes and for twin I and and for twin II are 
also shown. Solid lines represent contributions from twins I and II 
resulting in dashed line. 


susceptibility with wide maximum around w 60 K. 

Temperature dependence of magnetic susceptibility 
anisotropy measured in magnetically ordered state below 
Tat Ri 8.5 K is shown in Fig. 4. Angular dependence of 
torque multiplied by 2M,„oi/{niH^) [see Eq. (2)] measured 
at T = 4.2 K in the ac, a*b and be plane is shown in Figs. 
5, 6 and 7, respectively. Maxima i.e. minima of shown 
torque curves represent values of susceptibility anisotropy 
in those planes. Magnetic susceptibility anisotropy shown 
in Fig. 4 was measured by applying magnetic field in the 
direction where angular dependencies shown in Figs. 5, 6 
and 7 have minima. Below TV = 8.5 K large susceptibility 
anisotropy develops in the planes containing the b axis, a*b 
and be plane, while in the ac plane susceptibility anisotropy 
remains comparable to its value above Tjv. Susceptibility 
anisotropies AXia* =Xb- Xa* and AXic = Xb - Xc amount 
to —2.2 • 10^^ emu/mol at T = 4.2 K, which is slightly 
larger than the value of (x) at 7)v, see Fig. 1. In collinear 
antiferromagnet susceptibility anisotropy at T = 0 amounts 
to X± ~ X|| ~ where x± is expected to be temperature 
independent below TX [17]. Since T)v = 8.5 K in CuSb206, 
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FIG. 6. Angular dependence of torque multiplied by 2M,„gi/(mH^) 
[see Eq. (2)] measured in different fields in the a*b plane. Dashed 
line represents Eq. (2) expected when no spin reorientation takes 
place. Solid lines represent results of simulation discussed in text. 
Error hars are shown only when they are larger than the symbol. 


measured susceptibility anisotropies in planes with b axis will 
become even larger when T < 4.2 K. Observed magnitude 
of susceptibility anisotropy in comparison with {%) at Tn 
siginifies that in our sample (or at least the great majority of 
it) b axis is the easy axis. 

Angular dependence of torque multiplied by 2M,„oi / (mH^) 
measured at T = 4.2 K in the ac plane shown in Fig. 5 can 
be described by Eq. (2), although the signal is very weak 
even in the largest applied field of // = 8 kOe. If we fit 
the data to Eq. (2) we obtain the zeros of the sine curves 
in directions ri [101 ]* and c — a, with Xc-a > X[ioi]*- In 
Eig. 5 positions of macroscopic magnetic axes, mi and ma 
for twin I and II are also marked. Since our ESR results 
suggest that twins I and II make up for majority of the 
sample, we fitted the observed torque in the ac plane to 
T = A/sin(29 —20^)+A//sin(29 —20^/), where A/ and A// 

are amplitudes and 0^ and angles of axes and mi^ of 
twins I and II, respectively. The results of the contributions 
for each twin are shown in Eig. 5 by solid lines, and sum of 
those contributions by dashed line. 

Angular dependence of torque multiplied by 
2Mmoi/{niH^) in the a*b and the be plane measured at 
T — 4.2 K, shown in Figs. 6 and 7 respectively. Amplitude 
of T ■ {2Mfnot)/[mH^) does not change with applied field, 
signifying linear dependence of magnetization on magnetic 


0 50 100 150 200 0 50 100 150 200 



Angle (deg) Angle (deg) 

FIG. 7. Angular dependence of torque multiplied by 2M^gi/ 

[see Eq. (2)] measured in different fields in the be plane. Dashed line 
represents Eq. (2) expected when no spin reorientation takes place. 
Solid lines represent results of simulation discussed in text. Error 
bars are shown only when they are larger than the symbol. 

field. In low magnetic fields, H <5 kOe, angular dependence 
in both planes follows Eq. (2). In higher fields, however, a 
deviation from this dependence is observed. The deviation is 
most easily observed as a shift of the angle of maximum and 
minimum of torque away from the angle it is situated on in 
lower field. The deviation grows with applied magnetic field 
and it is slightly larger in the be plane. The observed deviation 
of torque is a result of reorientation of spin axis from the easy 
axis direction and is typical for uniaxial antiferromagnets 
when applied magnetic field becomes comparable to the 
spin flop fleld Hsf [20, 21]. In our case the deviation is 
observed in both a*b and be plane, although it is stronger 
in the latter. In the next section, we study influence of the 
magnetic anisotropy on the torque curves of this system using 
phenomenological approach in which magnetic anisotropy 
is described by the symmetry-obeying magnetocrystalline 
anisotropy energy. Comparison with measured data will 
allow us to determine the spin axis direction in both zero fleld 
and for//||fc > Hsf- 

B. Phenomenological approach to magnetic anisotropy in 
CuShaOfi 

Effect of spin reorientation phenomena in uniaxial antifer¬ 
romagnets on the torque curves has been successfully studied 
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FIG. 8 . Magnetocrystalline anisotropy energy allowed by symmetry 
in CuSb 206 , Eq. (4), with A'J > A";. Direction of magnetic axis m\, 
mi and m 3 with respect to crystal axes for each twin is given in table 

I. 


using phenomenological expression for uniaxial anisotropy 
energy [20, 21]. This approach has often been limited to only 
one plane in which the spins are allowed to rotate [21]. In 
this work we use the three-dimensional anisotropy energy al¬ 
lowed by symmetry where we allow the spin axis to rotate 
away from the easy axis in the direction which minimizes the 
total energy in finite magnetic field. This rotation is described 
by the rotation of the susceptibility tensor. From the rotated 
susceptibility tensor we can calculate the magnetization and, 
using Eq. (1), the torque. This approach was used recently 
to describe the measured torque curves in antiferromagneti- 
cally ordered phases of tetragonal Bi 2 Cu 04 in terms of easy 
plane anisotropy [ 22 ], and cubic CuaTeOg in terms of cubic 
anisotropy [23]. 

Magnetic susceptibility tensor for system with monoclinic 
space group P 12 i /«1 written in {a*,b,c) coordinate system is 
given by [18] 


0 Xb 0 . 

\Xae 0 Xc / 


(3) 


The above given expression for susceptibility tensor is 
valid for twins I and II. Assuming all twins have common 
a and b axes, for twins III and IV space group P2i/nll is 
more appropriate and unique axis is the a axis, as already 
discussed in Sec. IIC. The coordinate system used in that 
case is {a,b*,c) and the only nonzero nondiagonal element of 
the susceptibility tensor is then Xbc- 

Magnetic susceptibility tensor is diagonal in coordi¬ 
nate system spanned by the macroscopic magnetic axes 
{mi,mi,ms), where nii, mi and nis for each twin are given in 
table I. In what follows we use coordinate system spanned by 
magnetic axes {mi,mi,ms). 

The most general second-order magnetocrystalline 
anisotropy energy invariant to symmetry operations of 
monoclinic p phase of CuSb 206 can be written as [24] 

= ^1 sin^ 0 + sin^ 0cos2(|), (4) 


where Ki and K[ are phenomenological anisotropy constants 
and 0 and (|) are polar and azimuthal angles, respectively. To 
describe the observed result that axis mi is the easy axis {b) 


in CuSb 206 we set K[ > Ki. This gives the anisotropy en¬ 
ergy a global minimum at angles (0,(|)) = ( 7 t/ 2 , 7 t/ 2 ) to which 
we assign the direction of mi. The shape of this anisotropy 
energy is shown in Fig. 8 . The saddle point represents a di¬ 
rection in which the spins will reorient when magnetic field 
H > Hgp is applied in the easy axis direction mi. By symme¬ 
try requirements other two axes corresponding to saddle point 
and global maximum can be directed anywhere in the plane 
spanned by global magnetic axes mi and m 3 . Magnetic sym¬ 
metry, however, prefers these directions to be in the directions 
of the magnetic axes. This leaves us with two choices: pre¬ 
ferred direction of spins when they flop (the saddle point) is 
along nil or m 3 , the directions of minimal and maximal com¬ 
ponents of magnetic susceptibility in the {mi,ms) plane, re¬ 
spectively. We performed calculations for both and our results 
favor m 3 in the direction of the saddle point, such as shown in 
Fig. 8 . 

Magnetocrystalline anisotropy energy, Eq. (4), gives us 
the direction of the easy axis, i.e. spin axis in zero mag¬ 
netic field. In experiment finite magnetic field H(^,\(/) = 
//o(cos^sin\|/,sin^sin\|/,cos\|/) is applied, where Hq is the 
magnitude of the field, \|/ is the polar angle the field makes 
with axes m 3 and ^ is the azimuthal angle that the projection 
of H on the (mi,m 2 ) plane makes with the axes mi. The Zee- 
man term, can be written as 

= (5) 

In our calculations we use measured components of magnetic 
susceptibility expressed in emu/mol and magnetic field H ex¬ 
pressed in Oe. This gives ]7z expressed in erg/mol, so our 
anisotropy constants in expression (4) for anisotropy energy 
are given in units erg/mol. In low magnetic field, H Hsf, 
the susceptibility tensor in Eq. (5) is diagonal 


\dia ' 


^Xtrii 

0 

VO 


0 0 \ 

Xini 0 I 

0 Xni} / 


( 6 ) 


As magnetic field increases the spin axis start to reorient from 
the easy axis direction to minimize the total energy. Following 
procedure for in-plane reorientation in uniaxial antiferromag- 
nets [ 20 ] and more general reorientation in easy plane [ 22 ] and 
cubic [23] antiferromagnets, we describe the reorientation of 
spin axis as a rotation of the susceptibility tensor 


X{Q,^)=R{Q,^)-Xdia-'R^{QA), (7) 

where R(0,(|)) = R,,,, ((j) — il/2) • R,,,, (0 — Tl/2) is the matrix 
that describes the rotation of the susceptibility tensor [25]. To 
obtain the direction of the spin axis in finite magnetic field 
we minimize the total energy 

^total = ^a{QA) + ( 8 ) 

with respect to 0 and (|). This approach assumes collinear two 
sublattice antiferromagnet and it allows us to obtain the direc¬ 
tion of the spin axis as the direction of the minimal value of 
susceptibility tensor with respect to magnetic axes. 
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In our simulations we construct the diagonal susceptibil¬ 
ity tensor in such a way to reproduce the measured val¬ 
ues Xa, Xb and Xc in tensor (3) found in literature. Using 
Xmi = 1.99 • 10^^ emu/mol, Xm 2 = 2-0 • 10^^ emu/mol and 
Xnii = 2.39 • 10^^ emu/mol gives Xa* = 2.16 • 10^^ emu/mol, 
Xb = 2.0 • 10^^ emu/mol, Xc = 2.22 • 10^^ emu/mol and Xac = 
=F2 • 10^^ emu/mol, where — and -|- are signs for twins I 
and II, respectively. We then search for correct values of 
anisotropy constants Ki and K[ by calculating the depen¬ 
dence of magnetization on magnetic field and by finding val¬ 
ues of K\ and K[ which reproduce the value of the spin flop 
field Hsf « 1.25 T reported in literature from magnetization 
[4, 10, 12] and ESR measurements [4]. Using the follow¬ 
ing values for anisotropy constants Ki =5-10^ erg/mol and 
K[ = 2.21 • 10^ erg/mol gives the dependence of magnetization 
on magnetic field shown in Fig. 9 where spin flop is observed 
for magnetic field of H = 1.25 T applied along the m 2 = b 
axis. These values of anisotropy constants are of the order 
of magnitude typical for uniaxial antiferromagnet with this 
magnitude of Hsf [26]. In Fig. 9 we also show calculated 
magnetization for magnetic field applied in the a*, c, mi and 
ms direction. We observe that for > Hsf magnetization 
obtains the value it has for H\\ms. Absolute value of our cal¬ 
culated magnetization compares very well with the measured 
values from Ref. 4. 

The orientation of the spin axis obtained from our cal¬ 
culations for > Hsf is in the direction of m 3 , which 

= (—ti,0,y) for twin I, where r| = and y = 0.7584 

expressed in {a*,b,c) coordinate system. This means that 
spins have components along both a and c axis. If we ex¬ 
press the observed direction m 3 in the monoclinic {a,b,c) 
system, we get for spin axis direction {msa,msb,in 3 c) = 
(—0.6519,0,0.7712), where msa, msb and msc is the projec¬ 
tion of the m 3 axis on a, b and c axis, respectively. This gives 
|m 3 fl/m 3 c| « 0.85, which represents the ratio of magnetic mo¬ 
ment along a and c axes in our model. Obtained result com¬ 
pares rather well to the ratio of projections of magnetic mo¬ 
ment ma/nic « 0.8 1 observed in elastic neutron diffraction ex¬ 
periment performed on single crystal in // = 6 T applied along 
the b axis of CuSb 206 [10]. The orientation of spin axis for 
twin I in zero field and for > Hsf is shown in Fig. 10. 
Results for other twins are given in table I. Our result for the 
direction of spin axis is of course a consequence of the choice 
of the shape of phenomenological magnetocrystalline energy 
(4). To further prove that this choice is correct, we continue 
to calculate the angular dependence of torque using the same 
anisotropy energy. 

Magnetization M = x • H, expressed in emu Oe/mol is 
calculated for applied field H from the susceptibility tensor 
X obtained from the minimization of total energy ( 8 ), as de- 
cribed above. Magnetic torque is then calculated using expres¬ 
sion T = ml{2M„j„i) M X H, where m and are mass and 
molar mass of the sample, respectively. This gives torque in 
units of erg or dyn cm. We compare our measured torque with 
calculations for twin I only, since twins I and II are in great 
majority in our sample. Calculations for twin II for planes 
a*b and be in which the measurements were performed give 
the same torque curves as twin I due to symmetry, although 



FIG. 9. Dependence of magnetization on magnetic field obtained 
from simulations for twin I. 




FIG. 10. Orientation of spin axis in twin I in zero magnetic field and 
inH > Hsf applied along easy axis direction obtained from our phe¬ 
nomenological treatment of magnetic anisotropy in CuSb 206 . Spin 
axis are represented as double-headed arrows since the direction of 
spins at specific sites and with respect to one another cannot be de¬ 
termined by this method. For other twins the direction of spin axes 
with respect to crystal axes is given in table I. 


the direction of spin axis in the flopped phase is different (see 
table I). In Figs. 6 and 7 the results of calculations plotted by 
solid line are compared to the measured values. We multiply 
both measured and calculated torque by {2Mi^,i)/(inH^). The 
agreement of simulation and experiment is very good. Angu¬ 
lar dependence given by expression Fq. (2) (dashed line in 
Figs. 6 and 7), expected when no reorientation of spin axis 
takes place, describes the data well only in fields H <6 kOe. 
Having assumed that axis mj is in the direction of the sad¬ 
dle point of anisotropy energy, we obtain the result in which 
the deformation of the angular dependence in higher fields is 
slightly larger in the be plane than in the a*b plane. Calcula¬ 
tion with mi axis in the direction of the saddle point gives the 
opposite result which is not in agreement with the experiment. 
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IV. DISCUSSION 

Our study of magnetic anisotropy of magnetically ordered 
state of CuSb 206 was motivated by different reports on 
the direction of the easy axis in this system [5, 8-10, 12]. 
Torque magnetometry is sensitive to magnetic susceptibility 
anisotropy, see Eq. (2), and as such is a reliant experimen¬ 
tal method for determining the direction of easy axis in sim¬ 
ple collinear antiferromagnets [20, 21]. We characterized our 
sample by performing ESR measurements to determine the 
presence and type of crystallographic twins which can have 
two different directions of the elongation in two crystallo- 
graphically equivalent CuOg octahedra, as summarized in ta¬ 
ble I. Twins of type I and II with b axis as unique axis were 
shown to be dominant, and twins III and IV are found to make 
up less than 10% of the sample. The dominance of twins I and 
II was also reported in literature [3, 4, 14]. Erom our torque 
magnetometry measurements we determined that in our sam¬ 
ple, for which majority of twins are of type I and II, easy axis 
is along b axis. 

Direction of easy axis in twins III and IV depends on 
weather these twins share common crystal a and b axes with 
twins I and II. There is a possibility that for twins III and IV a 
axis is in the direction of the b axis of twins I and II. Easy axis 
of twins III and IV are then also in the b direction, but will 
be observed when field is applied in the direction of a axis in 
twins I and II. This scenario trivially explains how easy axis 
behavior can be observed for field applied along both a and 
b axis. However, there is another explanation we put forward 
in our analysis which assumes that a and b axes are common 
for all twins, but unique crystallographic axes of monoclinic 
crystal structure are different. In twins I and II b axis is unique 
axis of monoclinic crystal structure, and in twins III and IV it 
is a axis. Twins I and II thus have P12i/«1 symmetry while 
twins III and IV f’2 1 /«11 symmetry. This conclusion is drawn 
from the analysis of and x tensors in relation to the Neu¬ 
mann’s principle. We proceed with discussion assuming this 
scenario is realized in CuSb 206 . 

Magnetocrystalline anisotropy of twins III and IV must 
have the same shape as for twins I and II shown in Eig. 8 
when expressed in their respective macroscopic magnetic 
axes. Since it must be invariant to symmetry operations of 
the space group in the same way as magnetic susceptibility 
tensor, the 1112 axis in Eig. 8 must correspond to the unique 
axes of monoclinic crystal structure, i.e. m 2 = d. mi and hit, 
are then defined in the same way as for twins I and II, but are 
constrained to the be plane. Consequently, for twins III and 
IV easy axis is the a axis and in > H^f the spins flop in 
the direction m 2 , which is now in the be plane (see table I). In 
our sample twins III and IV together amount to less than 10% 
of the total sample, so we cannot quantitatively confirm this 
result. However, it would explain why in some cases reported 
in literature significant decrease of is observed along with 
the decrease of Xfo- Directions of easy axis and spin axis in 
flopped phase for each twin obtained from our phenomeno¬ 
logical model are summarized in table I. 

Eor one sample of CuSb206 it was recently reported that in 
the AEM ordered state easy axis and spin flop was observed 


for both a and b directions [12]. Indeed, the reported de¬ 
crease of susceptibility along both a and b axis in that sam¬ 
ple was comparable in magnitude corroborating this claim. 
Our analysis would explain this observation if different twins 
were present in this sample, regardless of the reason (com¬ 
mon or rotated a and b axes) but authors claim that they 
managed to completely detwin their sample. In that case 
the anisotropy of their sample should be described by dif¬ 
ferent type of anisotropy energy, such as found in tetragonal 
easy plane antiferromagnet Bi 2 Cu 04 [22]. The symmetry of 
this anisotropy energy is higher than required for monoclinic 
CuSb206. In more recent papers by the same authors on, pre¬ 
sumably, same samples the decrease of Xa was significantly 
smaller than of Xfe [13] and some twinning of the samples 
was reported [14]. The polishing of the surface which was 
used to detwin the samples in Ref. 12 and the method of tak¬ 
ing the Laue images to determine the presence of twins was 
shown to be insufficient by Prokofiev et al. who also used 
this procedure and then went on to further investigate the ap¬ 
parently detwinned samples by X-ray counter diffractometer 
only to show that twinning still persists [3]. In the same pa¬ 
per, Ref. 3, different attempts to completely detwin the sam¬ 
ples were reported to be unsuccessful and the conclusion was 
drawn that polishing only detwinns the surface. This adds to 
our claim that if two directions of easy axis, a and b, appear in 
the sample of CuSb 206 , they must come from different crys¬ 
tallographic twins. In one single crystal neutron diffraction 
study, which determined the b axis as easy axis, weak peaks 
from different twins were also observed, but these were, at 
least partially, removed from analysis. 

In order to unambiguously determine that the direction of 
easy axis can be different in different type of twin, as our anal¬ 
ysis suggests, a systematic study of the contribution of differ¬ 
ent twins should be performed in samples where both a and b 
axis were observed as easy axis. However, such quantification 
of the contribution of different twins was not reported so far 
in literature, except in this study and in Ref. 4, where samples 
with dominant twins I and II are found, and where easy axis is 
the b axis. A fully systematic study with different diffraction 
techniques, ESR and local magnetic techniques such as nu¬ 
clear magnetic resonance, should also include the systematic 
pressure treatment of the samples to determine what is the fac¬ 
tor that influences the contribution of different twins. While 
such analysis is beyond the scope of this work, it presents an 
interesting topic for future research. 

Einally, let us briefly comment on microscopic magnetic 
anisotropy that drives this system to specific arrangement of 
spins found in the ordered state. As stated in the introduc¬ 
tion, for Cu^+ spin 5=1/2 systems anisotropy of the mag¬ 
netically ordered states is driven by the weak anisotropic ex¬ 
change. In CuSb206 antisymmetric anisotropic exchange i.e. 
Dzyaloshinskii-Moriya interaction [27], is not allowed be¬ 
tween spins on chain and collinear arrangement of spins found 
in ground state suggests that it can be ignored. Symmetric 
anisotropic exchange then seems to be the reason behind the 
observed arrangement of spins. However, it was recently re¬ 
ported that spin-orbit coupling can induce weak single-ion 
anisotropy in spin 5=1/2 transition metal ions and that this 
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can be the driving mechanism behind arrangements of spins 
observed in some systems where standard approach fails [28]. 
Perhaps small spin flop field of 1.25 T observed for this sam¬ 
ple is connected to the gap induced by this type of weak 
anisotropy. 

V. CONCLUSION 

The antiferromagnetically ordered state of monoclinic 
CuSb 206 , which is still not satisfactorily understood, was 
studied using torque magnetometry as a sensitive probe of 
macroscopic magnetic anisotropy. Multiple crystallographic 
twins were detected by structural and ESR analysis which 
showed that in our sample dominant twins have PYlxjnX sym¬ 
metry with h axis as unique axis. Torque measurements reveal 
that h axis is the easy axis in this sample. However, our anal¬ 
ysis of crystal and magnetic symmetry of different twins in 
monoclinic CuSb 206 shows that both a and b axes can emerge 
as easy axis in different twins. Combining torque results with 
phenomenological approach to magnetocrystalline anisotropy 
energy we predict that when magnetic held H > Hsf is ap¬ 
plied along easy axis the spins flop in the direction of the 


macroscopic magnetic axis along which the susceptibility and 
tensor have maximal values and that this direction is dif¬ 
ferent for different twins. Our result nicely demonstrates how 
torque magnetometry can be used to determine the values of 
the spin flop field and the symmetry of the antiferromagneti¬ 
cally ordered state in both H < Hsf and H > Hsf by perform¬ 
ing measurements in fields which are significantly smaller 
than Hsf- Results presented in this work offer possibility to 
reconcile different reports on easy axis direction found in lit¬ 
erature and present a challenge for future studies of the crystal 
symmetry and magnetically ordered state of this system. 
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